Measuring storage and loss moduli using optical tweezers: broadband microrheology 
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We present an experimental procedure to perform broadband microrheological measurements 
with optical tweezers. A generalised Langevin equation is adopted to relate the time-dependent 
trajectory of a particle in an imposed flow to the frequency-dependent moduli of the complex fluid. 
This procedure allows us to measure the material linear viscoelastic properties across the widest 
frequency range achievable with optical tweezers. 
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In 1986 Ashkin and colleagues reported the first ob- 
servation of what is now commonly referred to as optical 
tweezers: a tightly focused beam of light capable of hold- 
ing microscopic particles stable in three dimensions 
Since then, several studies have adopted this technique as 
a tool for purposes as varied as trapping solid aerosols , 
measuring the viscosity of biomaterials ^, ,4] , the forces 
exerted by single motor proteins Q and the compliance 
of bacterial tails , or stretching single DNA molecules 
0. However, there remain a number of issues when opti- 
cal tweezers are used for microrheological measurements. 

Microrheology is a branch of rheology having the same 
principles as conventional bulk rheology (i.e. to study the 
linear viscoelastic behaviour of materials), but working 
on micron length scales. The linear viscoelastic proper- 
ties of a material can be represented by the frequency- 
dependent dynamic complex modulus G*(ljj), which pro- 
vides information on both the viscous and the elastic na- 
ture of the material. The conventional method of mea- 
suring G*{lo) is based on the imposition of an oscillatory 
stress a{uj,t) and the measurement of the resulting os- 
cillatory strain ^(io^t), or vice versa. The amplitudes of 
its in-phase and out-of-phase components are both pro- 
portional to the stress amplitude, with constants of pro- 
portionality defining, respectively, the storage (elastic) 
G'{u)) and the loss (viscous) G"{uj) moduli Q. 

Optical tweezers have been successfully used with New- 
tonian fluids for rheological purposes such as determining 
the fluid viscosity with high accuracy, measuring the hy- 
drodynamic interactions between particles or estimating 
the wall effect on the Stokes drag coefficient (i.e. Faxn's 
correction), as reviewed in Ref. 0]. Gonversely, when 
optical tweezers are adopted for measuring the viscoelas- 
tic properties of complex fluids the results are limited to 
the material high frequency response, discarding the es- 
sential information related to long times scales (i.e. low 
frequency) material behaviour. 
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The aim of this letter is to present a self-consistent pro- 
cedure for measuring the linear viscoelastic properties of 
materials, from non-oscillatory measurements, across the 
widest frequency range achievable with optical tweezers. 
In particular, the procedure consists of two steps: (I) 
measuring the thermal fluctuations of a trapped bead 
for a sufficiently long time; (II) measuring the transient 
bead displacement, from the optical trap centre, in re- 
sponse to a uniform fluid flow field entraining the bead. 
The flow is instantaneously switched on at time zero, 
by translating the whole fluid sample while the trap is 
held flxed. The imposed constant-velocity motion con- 
tinues until a steady displacement of the bead is reached. 
The analysis of the flrst step (I) provides: (a) the trap 
stiffness (k) - note that this has the added advantage of 
making the present method self-calibrated - and (b) the 
high frequency viscoelastic properties of the material, to 
high accuracy. The second step (II) has the potential to 
provide information about the viscoelastic properties of 
the material down to very low frequencies, limited only 
by the duration of the experiment. However, because 
of the harmonic nature of the optical trap, that tends 
not to trasmit high-frequency applied forces to the bead, 
the material's high-frequency response can not be deter- 
mined by this step. The full material viscoelastic spec- 
trum is thus resolved by combining the results obtained 
from steps (I) and (II). 

The new experimental procedure is analytically de- 
scribed through the analysis of the motion of a bead 
trapped in a stationary harmonic potential of force- 
constant K, where a uniform fluid flow field of magnitude 
Va can be exerted at time i = 0. The equation describ- 
ing the bead position r{t) V t can be derived by means 
of the generalized Langevin equation, which in three di- 
mensions is: 

md{t) = fnit)- f ((t-T){v{T)~Va{T))dT-Kr{t), (1) 
Jo 

where m is the mass of the particle, a{t) is its acceler- 
ation, v{t) is the bead velocity, Vs{t) is the fluid flow 
field velocity and /^(t) is the resultant of the stochastic 
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thermal forces acting on the particle. The integral term 
represents the viscous damping of the fluid, which incor- 
porates a generalized time-dependent memory function 

We now show how Eq. ^ evolves in the two cases 
mentioned above: when Vs{t) — and Vs{t) ^ 0, respec- 
tively. 

In the first case, where Vs{t) = 0, the optical tweezers 
can be calibrated by using the Principle of Equipartition 
of Energy: 
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(2) 



where, ks is Boltzmann's constant, T is absolute tem- 
perature and (r^) is the time- independent variance of the 
particle displacement from the trap centre, the origin of r. 
Despite all the possible methods for determining the op- 
tical trap stiffness (e.g . using the power spectrum or the 
drag force [lol. [ill. Il^) . the Equipartition method is the 
only method independent of the viscoelastic properties 
of the material under investigation and is thus essential 
from a rheological point of view. 

The thermal fluctuations of the trapped bead can also 
be investigated to determine the high frequency vis- 
coelastic properties of the material through analysis of 
the time dependence of the mean-square displacement 
(Ar2(r)) (MSD): 



(3) 



where t is absolute time and r is the lag-time. The av- 
erage is taken over all initial times t and the number of 
particles considered in the experiment, if more than one. 
In particular, using the assumptions adopted by Mason 
and Weitz in the study of the motion of thermally ex- 
cited free particles jl3| . at thermal equilibrium, where 

= and m {v{t)v{t)) = GksT V t, Equation 

([I]) yields, in the Laplace form, the velocity autocorrela- 
tion function: 
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ms + C(s) + k/s 



^.^(Af^Gs)), (4) 



where s is the Laplace frequency. Following Mason and 
Weitz [l^ in assuming that the bulk Laplace-frcqucncy- 
dependent viscosity of the fluid fi{s) is proportional to 
the microscopic memory function ({s) = 6TTafj{s), where 
a is the bead radius, Eq. Q can be written as: 



(Af2(s)) 



K 

ms 



(5) 



where the first term in the brackets reflects the viscoelas- 
ticity of the medium, the second term is related to the 
inertia of the bead and the third term takes into account 
the optical trap strength. It is easy to demonstrate that, 
for a micro-bead of density of order of Ig / cm? suspended 



in water, the product ms is negligible compared with the 
first term for the majority of the experimentally accessi- 
ble frequencies (i.e. s « 10^ s^^). With regard to the 
optical trap strength, two limiting cases can be distin- 
guished: (i) in the limit k/s — > 0, which can be obtained 
either for vanishing trap strength or for measurements 
performed at high frequencies, but lower than 10^ s""'^, 
Eq. ([5]) recovers the generalized Stokes-Einstein relation- 
ship derived by Mason and Weitz [l^l; (ii) in the limit 
k/s — > cx), which can be obtained either for a strong 
optical trap or for measurements performed at very low 
frequencies, Eq. (O gives the same result as if the bead 
were embedded in a purely elastic continuum with elas- 
tic constant of n/Q'Ka. For all intermediate cases, where 
0<k/s<oo, itis easy to show that, by analytic con- 
tinuation from Eq. ([5]), the complex modulus can be ex- 
pressed directly in terms of the time-dependent MSD: 
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(6) 



where ^Ar2(a;)y is the Fourier transform of (Ar2(r)). 

The second step of the procedure, which experimen- 
tally follows the first, consists of the analysis of the in- 
duced bead displacement from the trap centre due to an 
imposed time-dependent uniform fluid flow field Vs (t) en- 
training the bead. In this case. Equation ^ yields, in 
the Laplace form, the mean velocity of the particle: 



{v{s)) 



~as)Vs{s) 



ms + C(s) + k/s 



(r(s)) 



(7) 



where the brackets (...) denotes the average over several 
independent measurements (but not averaged over abso- 
lute time, since time-translation invariance has been bro- 
ken by the flow start-up aX t — Q). It is straight-forward 
to show that, by analytic continuation from Eq. ([7]), the 
complex modulus can be expressed directly in terms of 
both the imposed flow field and of the induced bead dis- 
placement from the trap centre: 



G*(c.) = sry(s)|,^. 



[k — muj^)iuj {fioj)) 



67ra 



{Vs{uj)-iu (f(w)} 



(8) 



where Vs{uj) and {f{ui)) are the Fourier transforms of 
Vs{t) and {r{t)), respectively. Note that, so far, the tem- 
poral form of Vs{t) is still undefined. Thus Eq. ([8]) repre- 
sents the general solution for G*{uj) independently of the 
temporal form of Vs (t) (e.g. sinusoidal function Vs sin(tt'i) 
or, as in this work, Heaviside step function VsH{t), where 
V^(t) = V t < and Vs{t) = 14 V i > 0). 

In principle. Equations ([6]) and ([8]) are two simple ex- 
pressions relating the material complex modulus G* (uj) 
to the observed time-dependent bead trajectory f{t) via 
the Fourier transform of either the r{t) itself (in Eq. 



3 



or the related MSD (in Eq. ([6])). In practice, the evalua- 
tion of these Fourier transforms, given only a finite set of 
data points over a finite time domain, is non-trivial since 
interpolation and extrapolation from those data can yield 
serious artifacts if handled carelessly. 

In order to express the two Fourier transforms in 
Equations ([6]) and ^ in terms of the N experimental 
data points (ifc, {Ar'^(T))k) and {tk, {r{t))^.), respectively, 
where fc = 1 . . . A^, which extend over a finite range, ex- 
ist only for positive t and need not be equally spaced, 
we adopt the analytical method introduced in Ref. p^ . 
In particular, we refer to Eq. (10) of Ref. [13] which is 
equally applicable to find the Fourier transform g{uj) of 
any time-dependent quantity g{t) sampled at a finite set 
of data points {tk,gk), giving: 



■uj 5 (w) = iujg{0) 

N 



E 

fc=2 



fffc - gfc-1 
tk — tk-l 



(9) 
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FIG. 1: (Color online) The normalised MSD vs. lag-time 
of a diameter bead in water and in two water based 

solutions of PAM at concentrations of 0.5 % w/w and 1 % 
w/w. The inset shows the same data, unnormalised. The line 
in the inset is the Einstein prediction of the MSD for a 
diameter bead in water at 25° C. 



where ^oo is the gradient of g{t) extrapolated to infi- 
nite time. Also g{0) is the value of g{t) extrapolated 
to t = 0+. Identical formulas can be written for both 
(^Ar^{uj)'^ and {f{uj)), with g replaced by (Ar^) and (f), 

respectively. This analytical procedure has the advan- 
tage of removing the need for Laplace /inverse-Laplace 
transformations of experimental data [l5| . 

We have tested Equations © and (O, via Eq. ©, by 
measuring both the viscosity of water and the viscoelas- 
tic properties of water-based solutions of polyacrylamide 
(PAM, flexible polyelectrolytes, M„ = 5 — 6x10'' g/niol, 
Polysciences Inc.) using optical tweezers as described be- 
low. 

Trapping is achieved using a CW Ti: sapphire laser sys- 
tem (M Squared, SolsTiS) which provides up to IWatt 
at 830nm. The tweezers are based around an inverted 
microscope, where the same objective lens, 100 x 1.3NA, 
(Zeiss, Plan-Neofluor) is used both to focus the trapping 
beam and to image the resulting motion of the parti- 
cles. Samples are mounted in a motorized microscope 
stage (ASI, MS-2000). Two CMOS cameras are used to 
view the sample, with bright-field illumination; one pro- 
vides a wide field of view (Prosilica EC1280M), while the 
other takes high speed images of a reduced field of view 
(Prosilica GV640M). These images are processed in real 
time at 2kHz using our own Lab VIEW (National Instru- 
ments) particle tracking software running on a standard 
desktop PC [13. 

The Brownian fluctuations of an optically trapped 
bead give rise to the time dependent (Ar'^(T)) shown 
in the inset of Fig. [TJ In the case of a bead immersed 
in a Newtonian fluid, it is expected that at short time 
intervals (thus small distances) the bead behaves as if it 
were free to diffuse. Indeed, the agreement between the 
observed (Ar^(T)) at short times of a trapped bead in 
water (circles) and the Einstein prediction for a freely 



diffusing bead (solid line) is good. As the time inter- 
vals increase the bead becomes influenced by the optical 
potential. This results in a plateau at large time inter- 
vals, where the (Ar^(T)) tends to 2(r^). It is interesting 
to note that the ratio of these two quantities (the MSD 
and twice the variance of the positional distribution) is 
a dimensionless parameter, independent of both the op- 
tical trap stiffness and the bead radius. It thus allows an 
explicit comparison between the dynamics of the fluids 
under investigation (Fig. [l}. Moreover, the onset point 
of the plateau region in Fig. [1] indicates the bottom limit 
of the frequency range within which the moduli can be 
determined by Eq. ([5]) , as for all the previous works using 
stationary optical tweezers. 

In Figure [5] we compare the responses of a 5fim diam- 
eter bead immersed in water (a Newtonian fluid) and in 
a water solution of PAM at 1% w/w (a non-Newtonian 
fluid), due to the imposition of a uniform fluid flow field 
having temporal behaviour as a Hcaviside step function 
Vs{t) — VsH{t), with different magnitude in the two mea- 
surements. Experimentally, the execution of a Heaviside 
step function is achieved by suddenly moving the mo- 
torised microscope stage at a predetermined speed and 
direction (here parallel to the x axis). The experiment 
runs until a steady displacement (Ax) of the bead from 
the trap centre is reached (i.e. until all the material's 
characteristic relaxation times are exceeded). In Figure^] 
the X component of the bead displacement has been nor- 
malised by Ax for a better comparison between the vis- 
coelastic character of the two samples. It is clear that 
while the Newtonian fluid reaches a steady value of the 
displacement almost instantaneously (as expected), the 
non-Newtonian fluid shows complex dynamics represen- 
tative of its viscoelastic nature. It is important at this 
point to note that, because of the harmonic nature of 
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the optical potential, at early times (i.e. for i — > or 
equivalently for uj — > oo), the trapping force exerted on 
the bead is actually small (i.e. Kf{t) — * 0) and the parti- 
cle moves almost at the same speed as the imposed flow 
(i.e. v{t) = Vs{t))] this implies that Eq. IH) becomes 
undefined at high frequencies. 
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FIG. 2: (Color online) The coordinates of a 5/im diameter 
bead vs. time for two different solutions and for two uniform 
fluid flow fields of difi'erent magnitude Vs at 25° C. In both 
the cases the data were averaged over three measurements and 
the X coordinate has been normalised by the steady state dis- 
placement Aa;. In water: Vs ~ 20 fim/s and Aa; — 0.523/^m. 
In 1% w/w of PAM, Vs = 3fim/s and Aa; = f .f55^m. The 
inset highlights the start-up behaviour of both the above sys- 
tems. 

The broadband microrheological measurement with 
optical tweezers is achieved by combining the frequency 
responses obtained from both the methods introduced 
above. In particular, the material's high frequency re- 
sponse is determined by applying Eq. © {via Eq. © 
with (Ar^)fc replacing gk) to the (Ar^(r)) measurements; 
whereas, the low frequency response is resolved by apply- 
ing Equation ([8]) {via Eq. ^ with replacing gk) to 
the data describing the bead's transient response to the 
motion of the stage. 

A typical result of this procedure for a non-Newtonian 
fluid is shown in Figure [3] while, in the case of water, a 
constant viscosity of ry = 8.69 x 10^"* ± 6 x lO^^Pa ■ s is 
measured over five frequency decades at 25°C. It is evi- 
dent that, although there is some noise in the frequency 
domain that has propagated from genuine experimental 



noise in the time-domain data, there is a clear overlap- 
ping region of agreement between the two methods that 
makes the whole procedure self-consistent. Moreover, it 
confirms the ease with which the low-frequency material 
response can be explored, right down to the terminal re- 
gion (where G' oc uj'^ and G" oc w), which is the current 
limitation for microrheological measurements performed 
with optical tweezers. However, in order to remove the 
genuine noise, a simple smoothing operation of the orig- 
inal data is sufRcient and the results are shown in the 
inset of Figure El 
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FIG. 3: (Color online) Storage (squares) and loss (circles) 
moduli vs. frequency of a solution of 1% w/w of PAM in 
water measured by means of both Eq. ((6} (solid symbols) and 
Eq. (|8]) (open symbols) applied directly to the experimental 
data presented in Fig. [T]and Fig. [2l respectively. The inset 
shows the moduli of the same solution as above, with both 
Eq. ([6} and Eq. ([8)l applied to the original data mentioned 
above, but smoothed. 

In summary, we have presented a self-consistent and 
simple experimental procedure, coupled with an analyt- 
ical data analysis method, for determining the broad- 
band viscoelastic properties of complex fluids with opti- 
cal tweezers. This method extends the range of the fre- 
quency response achieved by conventional optical tweez- 
ers measurements down to the material's terminal region. 
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